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1 Introduction 

High energy electron-positron or Bhabha scattering [1] is among the most important and 
carefully studied processes in particle physics. It provides a very efficient tool for luminosity 
determination at electron-positron colliders and thus mediates the process of extracting 
physical information from the raw experimental data [2-4] . In particular it is crucial for the 
high precision physics program at a future International Linear Collider (ILC). The small- 
angle Bhabha scattering is particularly effective as a luminosity monitor while the large- 
angle scattering is suggested as a tool to disentangle the luminosity spectrum at the ILC [5, 
6] . Bhabha scattering involves stable charged leptons both in the initial and the final states 
and, therefore, it can be measured experimentally with very high precision. On the theory 
side it can be reliably computed order-by-order in perturbation theory. These properties 
make Bhabha scattering an ideal "standard candle" for electron-positron colliders. Since 
the theoretical accuracy directly affects the luminosity determination, remarkable efforts 
have been devoted to the study of the radiative corrections. The one- loop corrections have 
been known in the full electroweak theory for several decades [7]. Recently the calculation 
of the two- loop QED corrections has been completed [9-20]. This result is sufficient for 
the small-angle or low-energy process. For the large-angle scattering above the electroweak 
scale the electroweak corrections become equally important and must be taken into account 
to match the designed experimental accuracy of a few permill [8]. The full evaluation of 
the electroweak two-loop corrections remains beyond the reach of available perturbative 
techniques. However, once the center-of-mass energy y/s is far larger than the electroweak 
gauge boson mass M^z, the cross section receives virtual corrections enhanced by powers of 
"Sudakov" logarithm ln(s / Myy z ) which dominate the two-loop corrections. Over the last 
few years the electroweak Sudakov logarithms have been studied within different approaches 
and the corresponding high-order corrections to various processes have been obtained [21— 



34]. 
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In this paper we generalize the next-to-next-to-next-to- leading (N 3 LL) result of Ref. [30] 
for the // —¥ f'f four-fermion annihilation process to Bhabha scattering. The result 
includes all the two-loop logarithmically enhanced corrections. The quartic, cubic and 
quadratic logarithms are obtained in full SUl(2) x f/y(l) electroweak standard model, 
while the linear logarithmic term is computed in a simplified SUl{2) theory sufficient for 
the practical applications. In the next section we outline the approach and present the 
analytical one- and two- loop results for the differential cross section in an SUl{2) model. 
In Section 3 we present the numerical result for the cross section in full electroweak the- 
ory. The main difference in respect to the analysis given in Ref. [30] is in the presence of 
both annihilation and scattering contributions to the Bhabha cross section and in using 
the electron mass as a collinear regulator for QED radiative corrections. 

2 High-energy asymptotic of radiative corrections in SUl{%) model 

We study the electron-positron scattering in the high-energy and wide-angle limit when 
all the kinematical invariants are of the same order and far larger than the gauge boson 
mass, |s| ~ |t| ~ \u\ 3> M% w . In this limit the asymptotic energy dependence of the field 
amplitudes is dominated by Sudakov logarithms [35, 36] and governed by the evolution 
equations [37-40] . The method of the evolution equations in the context of the electroweak 
corrections is described in detail for the fermion pair production in Ref. [30]. The evolution 
equations for the scattering (annihilation) amplitude are linear which results in the expo- 
nentiation of the Sudakov logarithms. To outline the general structure of the corrections 
let us for a moment neglect the hypercharge interaction and consider a spontaneously bro- 
ken SUl{2) gauge model with the gauge boson of mass M and twelve massless left-handed 
fermion doublets. This model retains the main features of the massive gauge boson sector 
of the Standard Model. In this case the result can be presented in a simple analytical 
form and constitutes the basis for the further extension to the full electroweak theory. The 
amplitude can be factorized into the following form 

^JSW^Vl, (2.1) 

where g([i 2 ) is the SUl(2) coupling renormalized in MS scheme, Q is the Euclidean mo- 
mentum transfer Q 2 = —s for annihilation and Q 2 = —t for the scattering process, J- is 
the electron vector form factor, and the reduced amplitude A. is a vector in an isospin and 
chiral basis. The solution of the evolution equations read 



T = F (a(M 2 ))exp 




M 2 X 



a ^7(«(^ , )) + c(«(^)) + e(«(M 2 )) 

M 2 x 



(2.2) 



A = Pexp 



M 2 X 



Ao(a(M 2 )) , (2.3) 



where the perturbative functions Fo(a), 7(a) etc. are given by a series in the coupling 
constant a(/U 2 ) = g 2 {[i 2 ) / (4?r) . By calculating the functions entering the evolution equa- 
tions order by order in a one gets the logarithmic approximations for the amplitude. In 
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the SUl{2) model with the degenerate Higgs and gauge boson mass the parameters of the 
Eqs. (2.2,2.3) are known to the N 3 LL approximation which includes all the two-loop loga- 
rithmic corrections [30] . By analytical continuation and crossing symmetry it is straightfor- 
ward to find the electron-positron scattering and annihilation amplitudes for the physical 
values of the Mandelstam variables and the corresponding contributions to the cross sec- 
tion. 



2.1 Differential cross section 

The perturbative series for the Bhabha cross section is of the following form 
da a 2 (s)(l-x) 2 



dfl 64s 
- 2x ( 1 + 



x- 

a 
4tt 



4tT b 167T 2 



a d4 2) 



+ !l + 4>" )+ 16. 2 



° dai 2) 



— da {2) 



(2.4) 



where the variable x = (1 — cos 9)/ 2 is a function of the scattering angle 9, a = a(M 2 ), and 
the annihilation S, scattering T, and interference ST contributions are given separately. 
The coupling constant in the Born approximation is normalized at fj? = s and can be 
expanded as follows 



as 



a 



1 



a 



t /° ln (i) + 



a y 
4tt) 



(2.5) 



with (3q = 19/6 and j3\ = —35/6. The corrections can be further split according to the 
power of the large logarithms 

d^)=d^ln^(^) + d.«ln(^) + d^, 



da^ 



d ^' ln \W 2 ) +d ^ 



\M 2 ) 



daf } In 



M 2 ) 



Oil) 



(2.6) 



where derp are functions of x, da^ corresponds to the one-loop leading logarithmic (LL) 
correction, da£ corresponds to the one-loop next-to-leading logarithmic (NLL) correction, 
and so on. The high-energy asymptotic of the one-loop corrections reads 



da 
da 

da 



(i) 

S 2 

(1) 
S 1 

(1) 
S 



9 + 21n(l 

235 
1~8~ 



101n(x) 



2 51n(g) 

6TT + 

1 — X 



5(1 - 2x)ln 2 (x) 



2(1 



(2.7) 



da. 
da. 



da, 



(i) 

T 2 

(1) 
T 1 

(1) 
T 



-3 , 

9 + 21n(l -x) + 21n(x) , 
935 

+ 2vr 2 + 21n(l-x)ln(x) 



18 

(8 + 7x) ln(x) + 5(2 - 2x + x 2 ) ln 2 (x) 



3(1 



2(1 -x) 2 



(2.f 
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da 



(1) _ _o 
ST 2 — ° ' 

(1) 



do£f, ! = 9 + 21n(l 



41n(x) 



da 



(i) 

ST 



1 



7T 2 + ln(l - x) ln(x) 



235 

~18~ 2 

(23 + 7x)ln(x) 5(l + x 2 )ln 2 (a;) 



+ 



6(1 - x) 



4(1 - x) 2 



(2.9) 



By reexpanding the exponents (2.2, 2.3) one gets the two-loop logarithmic corrections which 
are given by Eqs. (A.1-A.3) in the Appendix A. Eqs. (2.7, A.l) also give the corrections to 
the differential cross section of the e + e~ — > [i + process or, in general, of a different flavor 
same-isospin fermion pair production. For the production of an opposite-isospin fermion 
pair such as e + e~ — > vv the one- loop coefficients read 



da 
da 

da 



(i) 

S 2 

(1) 
S 1 

(1) 
S 



9 - 101n(l - x) + 21n(x) 



235 2 

3tt 2 

18 



ln(x) _!_(!- 2x)ln\x) 



l-x ' 2(1 -xf 
and the two-loop coefficients can be found in Appendix A. 



(2.10) 



2.2 Total annihilation cross section 



The differential cross section of the pure annihilation process can be integrated over solid 
angle to get the total cross section a = ^ n (a/47r) n o"( n ), where a^ = 7ra 2 (s)/(48s) is the 
total cross section in the Born approximation. The one- and two-loop corrections to the 
total cross section are 



a 



80 



3 lU (w) 



a 



(2) 



Mi) 



461 



„ >" 3 (i) 



/85 2 

U + 3 ' 

/5875 



a 

37ir 2 



(0) 



>" 2 (i) 



+ 



16019 



for e + e 



216 
and 



^^7T 2 - 122C(3) + 15\/37r + 26^012 / " 
18 



a' 



Pill) 



a 



(i) 



a 



(2) 



"Mi) 

Mi)- 

/ 24481 341 
+ V~21~6 18" 



26, / s \ 
137, 3 / s \ 



3tT< 



158 

~9~ 
451 37vr 2 
~9 3~~ 



a 



(0) 



>° 2 (i) 



vr 2 - 122C(3) + 15\/3vr + 26\/3Cl 2 (|)^ In (-^) 



<7 



(°)(2.12) 



for e + e — >• z^i/ process. Here CI2 (§) = 1.014942 ... is the value of the Clausen function and 
C(3) = 1.20206. . . is the value of the Riemann zeta-function. Similar formulae have been 
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derived in Refs. [25, 29, 30] for a model with six massless left-handed fermion doublets. 1 
Our choice of the fermionic content provides more accurate model of the electroweak theory 
with three generations of fermions. 



3 Two-loop electroweak corrections to Bhabha scattering 

In the full electroweak standard model the gauge coupling depends on the electron chirality 
and the cross section can be written as follows 

^ = £ E (^f J + daf J + 2dafJ), (3.1) 

I,J£{L,R} 



where the index / (J) stands for the initial (final) electron chirality, the MS SUl(2) 
coupling constant a = a e /sin 2 Ow is a ratio of the QED coupling constant a e an the weak 
mixing angle 0\y normalized at fi = Miy. In the Born approximation one has 

d-f L (0) = (1 + %) 2 (1 " x) 2 , d^ (0) = ^(l - x)\ daf R (0) = d< (0) = %* 2 , 



A n T (0) _ (l , t 2 w \ 2 (1-z) 2 j T (0) _ ,4 (I-*) 2 j T (0) _ j T (0) _ t 4 w I (o 2) 
d(J 5T(0) _ _ A , tjy\ 2 (l-xf , ST (0) _ ,4 (1^)! dfJ 5T _ d(J ST _ n 

where t 2 ^ = tan 2 #iy. The total electroweak Born cross section of Bhabha scattering then 
reads 

j Born 2 , . 

= zr—zd 1 - *) C 1 + *w0 2 + 8t wtt - 8x + 12x2 - 8x3 + 3a;4 ) ) • ( 3 - 3 ) 

We replace the hypercharge coupling constant by a tan 2 9w and consider a single-parameter 
perturbative series for the cross section 



da da BoTn a 2 
+ 



dn dU Asx 2 



+ *?>^(4) + d^h.(^r) +0 (l)| +0 (O 



(3.4) 



The analysis of Sudakov logarithms in electroweak standard model with the spontaneously 
broken SUl{2) x £/y(l) gauge group is similar to the treatment of the pure SUl(2) case 
considered in the previous section. The main difference is in the presence of the massless 
photon. The photon contribution to the virtual corrections results in the QED Sudakov 
logarithms of the form lnQ 2 /A 2 and hiQ 2 /m 2 , where rrif is a light fermion mass and A is 
the auxiliary photon mass introduced to regulate the infrared divergences. By factorizing 



Eq. (3) [29] and Eq. (44) [30] in the coefficient of the single logarithmic term 48049/216 - 1679tt 2 /18 
should be replaced by 38221/216 - 1601tt 2 /18. In Eq. (4) [29] and Eq. (45) [30] in the coefficient of the 
single logarithmic terms 38005/216 - 383tt 2 /18 should be replaced by 28177/216 - 305tt 2 /18. 
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the QED logarithms from the solution of the full electroweak evolution equations one gets 
pure electroweak logarithms of the form ln(Q 2 /M^ z ). The QED logarithms in turn form 
an exponential factor. In Ref. [30] this factor has been evaluated to N 3 LL accuracy for 
massless fermions. For our analysis we need only the NNLL approximation. In the limit 
m 2 <C A 2 <C Myyz ^ s ' w here / 7^ t stands for all the fermions except the top quark, for 
the annihilation amplitude it takes the following form 



U 



exp 



x W 



q e (A 2 
2tt 



+ 13 + 2 In 

+ o( 



a 2 In 



76 8 , 

h - In 

27 9 



1 — x 

X 

s 

A 2 



In 



7T 

8 a f 
27 



No 



(3.5) 



Here N g = 3 stands for the number of generations and —1/2 subtracts the top quark con- 
tribution. A non-standard limit with mj C A has been used to facilitate the procedure of 
separating the QED and electroweak logarithms since it allows for the calculation with the 
unbroken theory fields. The normalization of the QED factor depends on the factorization 
scheme and may include nonlogarithmic corrections in ot e . We chose the normalization 
U(s = A 2 ) = 1. The corresponding expression for the scattering amplitude may be easily 
obtained by crossing symmetry and analytical continuation of Eq. (3.5). 

The infrared divergent virtual QED contribution should be combined with the real 
photon radiation to get a physically observable infrared finite cross section. In practice, 
the massive gauge bosons are supposed to be detected as separate particles. Thus it is of 
little physical sense to treat the hard photons with energies far larger than Mz,w separately 
because of gauge symmetry restoration. In particular, the radiation of the hard real photons 
is not of the Poisson type due to its non-Abelian SU(2) l component. Therefore, we restrict 
the analysis to semi-inclusive cross sections with the real emission only of photons with 
energies far smaller than MzWi which is of pure QED nature. The standard Monte Carlo 
event generators [41-43] adopt a different limit on the photon mass A < mj. Let us now 
consider the evolution of the QED factor (3.5) to the small values of A. At every threshold 
A = m,f the anomalous dimensions of the corresponding evolution equation change due to 
the decoupling of the massive virtual fermion and one has to match the two solutions which 
result in the logarithms of the form ln(s/m 2 ) in the exponent. At the lowest threshold 
A = m e there is an additional nonlogarithmic one-loop matching correction because for 
A <C m e the electron rather than photon mass regulate collinear divergence. This matching 
correction can be obtained by comparing the one-loop corrections to the vector form factor 
in the two limits [15]. The result for the QED factor, which is now compatible with the 
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Figure 1. The one- and two-loop electroweak corrections in % to the Born cross section as function 
of the scattering angle for (a) = 500 GeV and (b) y/s = 1 TcV. The shaded area corresponds 
to the region where Sudakov approximation is not reliable. 



Monte Carlo event generators, reads 



U 



1 + 



7T 



+3 In 



m- 



1 vr 2 

1 + T 



+ 2 In 



exp 



a e (m; 



In 



2tt 



A 2 



In 2 



4 +2 In JL 

m eJ V We 



1 In 



in: 



+ 
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7T 



771, 



/19 1, A - 
+ 18 + 3 ln — 



/ 



77T,j 



+ O [ai In 



(3.6) 



The sum of the virtual and real QED corrections to the cross section is finite for A = 
and gives a factor which depends on the Mandelstam variables, the fermion masses, and 
the experimental cuts, but not on Mz t w- 

Thus, in comparison to the SUl(2) model the evolution equations describing the elec- 
troweak Sudakov logarithms are modified by the presence of the hypercharge interaction 
and by the subtraction of the QED logarithms. All the parameters of the evolution equa- 
tions are known to the NNLL approximation in the full electroweak theory including the 
effect of the gauge boson mass splitting [25] and the corresponding one- and two-loop 
corrections can be easily obtained. We present the result in semi-numerical form for the 
following input values 

a~ 1 (M z )= 128, sin 2 9 W (M Z ) = 0.231, M z = 91.2 GeV , M w = 80.4 GeV , 

(3.7) 

and choose Mw as the argument of the Sudakov logarithm. The high-energy asymptotic 
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Figure 2. The one- and two-loop electroweak corrections in % to the Born cross section as 
function of the center-of-mass energy for (a) 9 — 50° and (b) = 90°. 



of the one-loop corrections read 

d4 1} = - 0.34 + 1.10a: - 1.66x 2 + l.lOx 3 - 0.34x 4 , (3.8) 

da{ 1] = 1.21 - 3.33x + 4.99x 2 - 3.33x 3 + 1.21x 4 

+ (0.10 - 0.52x + 0.78x 2 - 0.52x 3 + 0.10x 4 ) ln(l - x) 

- (0.02 - 0.49x + 2.05x 2 - 2.24x 3 + 0.75x 4 ) ln(x) , (3.9) 

dtrft = - 0.49 + 1.78a; + 2.60x 2 - 5.34x 3 + 2.61x 4 
+ 0.03 (x - x 2 + x 3 ) ln 2 (l - x) 

- (0.05 - 0.16x + 0.28x 2 - 0.16x 3 + 0.05x 4 ) ln(l - x) 
+ (0.31 - 0.24x - 0.60x 2 + 0.48x 3 ) ln 2 (x) 

+ (1.21 - 2.51x + 2.64x 2 - 0.97x 3 + 0.05x 4 ) ln(x) 

+ (0.10 - 0.44x + 0.42x 2 - 0.13x 3 ) ln(l - x) \n(x) . (3.10) 

The two- loop LL, NLL, and NNLL electroweak corrections are given by Eqs. (B.l) 
in the Appendix B. The N 3 LL electroweak correction is sensitive to fine details of the 
gauge boson mass generation and currently the exact result is not available. We evaluate 
them approximately within a pure SUl(2) model considered in the previous section. The 
numerical result is given by Eq. (B.2) in the Appendix B. Note that Eq. (B.2) correspond 
to the decomposition (3.4) where all the coupling constants are normalized at the gauge 
boson mass (c/. Eq. 2.4). In the full standard model the hypercharge interaction and gauge 
boson mixing result in a linear-logarithmic contribution which is not accounted for in this 
approximation. It is, however, suppressed by the small factor sin 2 6w ~ 0.2. Therefore, 
the approximation gives an estimate of the coefficient in front of the linear electroweak 
logarithm with 20% accuracy. As we will see, such an uncertainty in the coefficient of the 
two-loop linear logarithm is negligible for practical applications. 
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Figure 3. Separate two-loop logarithmic contributions in % to the Born cross section as function 
of (a) the scattering angle for s/s = 1 TeV and (b) the center-of-mass energy for 8 = 50°. The 
shaded area corresponds to the region where Sudakov approximation is not relaible. 



Let us now consider the numerical effect of the corrections. We plot the one- and two- 
loop electroweak corrections relative to the Born cross section for the phenomenologically 
interesting intervals of the center-of-mass energy and the scattering angle in Fig. 1 and 2. 
For the small and very large scattering angles the Sudakov approximation is not reliable 
since there either the t or u invariant is comparable with M^ vz an d the power suppressed 
terms become important. In Fig. 1 the shaded regions denote the kinematic domain where 
either t or u are less than 10M§ and the Sudakov limit used in this calculation begins 
to falter. The one-loop correction strongly depends on the center-of-mass energy and 
scattering angle and is large and negative in the practically interesting region. Overall 
the two-loop correction is positive and enters at the 1 percent level everywhere in the 
angular distribution, with the largest corrections coming at 9 ~ 90°. The breakdown of 
the two-loop corrections into the individual logarithmic contributions is given in Fig. 3. 
The corrections exhibit typical sign-alternating behavior with significant cancellation of 
the individual logarithmic terms. For example, at \fs = 1 TeV and = 50° the numerical 
structure of the corrections read 




-9.45% 

One loop: - 3.20% { NLL: 7.35% Two loops: 0.67% { 

-1.09% 



LL: 0.62% 
NLL: -0.41% 
N 2 LL: 1.36% 
N 3 LL: -0.90% 

(3.11) 

Thus a twenty percent uncertainty in the coefficient of the two-loop linear logarithm would 
only induce an uncertainty on the order of one permill in the cross section. 

To conclude, we have computed the dominant two- loop electroweak corrections to 
high-energy wide-angle Bhabha scattering. The corrections can be as large as 10% in one 
loop and 1% in two loops. Our result completes the perturbative analysis of the Bhabha 
scattering necessary for the luminosity determination at the ILC. 
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A Two loop results in SUl(2) model 



Here we present the analytical result for the two-loop logarithmic corrections to the cross 
section of Bhabha scattering in the SUl(2) model with the gauge and the Higgs boson of 
the same mass M. The annihilation contribution reads 



9 

2 ' 

(2) 143 



da£> 3 = — - 61n(l - x) + 301n(x) , 

d<4 2) 2 = + ll7f2 + y ln ( X - x) + 2 ln 2 (l - x) - 8 ln(x) ln(l - x) 

(535 - 445x) m(x) (91 - 182x + 76x 2 ) ln 2 (x) 



6(1 — x) 2(1 — x) 2 ' 

/125 50vr 2 101n(x) 5(1 - 2x) ln 2 (x) \ , . 
/ 1075 - 1250x (74- 148x + 38x 2 W 2 \ , , , 

+ (- 6d-,) + - — 3<i-»)' ) lnW 

^ (631 - 806x) ln 2 (x) 19(1 - 2x) ln 3 (x) 
4(1 -x) 2 + (1 -x) 2 ' 



(A.l) 



The scattering contribution reads 



da, 



(2) 



9 



'T4- 2 > 

143 

ddyl = 61n(l — x) — 61n(x) , 



d4 2) 2 = - ^ + 8vr 2 + — ln(l -x) + 2 ln 2 (l - x) - 14 ln(x) ln(l 



+ 



4 6 
(41 - 131x) ln(x) (26 - 22x + llx 2 ) ln 2 (x) 



6(1 -x) 2(1 -x) 2 

d 4 2) i = ^ " 122 C(3) + W3C1 2 (£) + 15V3vr + ^±M^ 
216 V 3 / 1 — x 



/125 8^ 2 (17-7x)ln(x) (10 - 30x + 15x 2 ) ln 2 (x) \ , . 

+ {-T + lT + i-x (i-x) 2 

. N1 . /1030- 505x (80 - 88x + 44x 2 )vr 2 \ , . . 

+ 41n(x) ln 2 (l - x) + ( lg(1 _ g) + 3(1 _ g)2 J M*) 

(414 - 502x + 263x 2 ) ln 2 (x) (10 - 22x + llx 2 ) ln 3 (x) 
+ 12(1 -x) 2 + (1-x) 2 ' 



(A.2) 
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The interference contribution reads 



da (2) - 9 

aa ST 4 — 2 ' 

143 

da^ 3 = — — - 61n(l - x) + 121n(x) , 

dag? 2 = - y ^ 2 + y ln(l - x) + 2 ln 2 (l - x) - 11 ln(x) ln(l - x) 

_ (247 - 157s) ln(x) _ (7 + 16x + 7x 2 ) ln 2 (x) 
6(1 -x) 4(1 -x) 2 ' 

dag? ! = ^ - 122 C (3) + 26^Cl 2 g) + 15^vr - 

125 llvr 2 (27-7x)ln(x) (27 - 64x + 27x 2 ) ln 2 (x) 



+ {-T + ^r + 2d -V 2(i-^ " )Mi--i 

, , o/ x /2195 - 3245x (13 + 41x)vr 2 \ , , , 



(49 + 32x - 431x 2 ) ln 2 (x) _ (8 + 2x 2 ) ln 3 (x) 
24(1 -x) 2 (1 -x) 2 



(A.3) 



For the e + e — >■ z^f annihilation cross section the two-loop logarithmic corrections read 



A ( 2 ) 9 
d 4 4 = 2' 

d4 2) 3 = — + 301n(l - x) - 61n(x) 



daf\ = + llvr 2 - ^ ln(l - x) + 38 ln 2 (l - x) - 8 ln(x) ln(l - x) 



^ (107 - 89x) ln(x) (1 - 2x + 4x 2 ) ln 2 (x) 



6(1 -x) 2(1 -x) 2 

d 4 V ^_ 122C( 3) + 2^0^ 5 ) +15 ^-2^! 

+ ( _ f + 3? _ 2 ^ + a^ )ln(i _ i) 

/215 - 250x (14 - 28x + 50x 2 )vr 2 \ , , , 

+ (^xT + 3(1 -x) 2 J HX) 

(11 -46x)ln 2 (x) (l-2x)ln 3 (x) , AJ . 

+ 77. \ 2 — + V jz 1 , 2 V 7 • A.4 

4(1 -x) 2 (1 - x) 2 
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B Two loop electroweak numerical results 

The two-loop electroweak logarithmic corrections to the Bhabha cross section to the NNLL 
approximation read 

daf ] = 0.34 - 1.21x + 1.81x 2 - 1.21x 3 + 0.34x 4 , 

daf ] = - 1.43 + 3.51x - 4.81a; 2 + 2.91x 3 - 1.07x 4 

- (0.16 - 0.81a; + 1.22x 2 - 0.81x 3 + 0.16x 4 ) ln(l - x) 

- (0.18 + 0.74x - 4.23x 2 + 4.91x 3 - 1.71x 4 ) ln(x) , 

daf ] = 5.78 + 6.42x - 20.04x 2 + 8.19x 3 + 3.98x 4 

+ (0.05 - 0.21x + 0.30x 2 - 0.21x 3 + 0.05x 4 ) ln 2 (l - x) 
+ (0.97 - 3.27x + 3.62x 2 - 1.48x 3 - 0.13x 4 ) ln(l - x) 

- (0.29 + 0.29x - 4.06x 2 + 5.22x 3 - 1.97x 4 ) ln 2 (x) 

- (3.74 - 7.56x + 9.03x 2 - 4.45x 3 + 0.43x 4 ) ln(x) 

- (0.88 - 2.76x + 2.13x 2 - O.lOx 3 - 0.28x 4 ) ln(l - x) ln(x) . (B.l) 

The numerical result for the two-loop N 3 LL corrections to the Bhabha cross section in the 
SUl(2) model reads 

daf } = - 0.30 + 0.87x + 0.41x 2 - 2.21x 3 + 1.24x 4 

+ (4.25 - 18.23x + 36.60x 2 - 35.50x 3 + 12.89x 4 ) ln(l - x) 

+ (0.63 - 0.38x + 1.88x 2 - 2.13x 3 ) ln 3 (x) 

+ (0.18 - 0.38x - 1.75x 2 - 0.03x 3 ) ln 2 (x) 

+ (18.97 - 10.32x + 6.39x 2 - 9.83x 3 - 5.21x 4 ) ln(x) 

+ (0.25 - 0.75x + 0.75x 2 - 0.25x 3 ) ln 2 (l - x) ln(x) 

- (0.63 - 3.56x + 5.25x 2 - 2.31x 3 ) ln(l - x) ln 2 (x) 

+ (0.27 - 1.60s + 2.40x 2 - 1.06x 3 ) ln(l - x) ln(x) . (B.2) 
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